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C^ ■ 

The differential migration of ions in an applied electric field is the basis for separation of 
chemical species by capillary electrophoresis. Axial diffusion of the concentration peak 
limits the separation efficiency. Electromigration dispersion is observed when the con- 
centration of sample ions is comparable to that of the background ions. Under such 
conditions, the local electrical conductivity is significantly altered in the sample zone 
making the electric field, and therefore, the ion migration velocity concentration de- 
pendent. The resulting nonlinear wave exhibits shock like features, and, under certain 
O ' simplifying assumptions, is described by Burgers' equation (S. Ghosal and Z. Chen Bull. 

'o , Math. Biol. 2010 72, pg. 2047). In this paper, we consider the more general situation 

where the walls of the separation channel may have a non-zero zeta potential and are 
therefore able to sustain an electro-osmotic bulk flow. The main result is a one dimen- 
sional nonlinear advection diffusion equation for the area averaged concentration. This 
homogenized equation accounts for the Taylor- Aris dispersion resulting from the vari- 
ation in the electro-osmotic slip velocity along the wall. It is shown that in a certain 
range of parameters, the electro-osmotic flow can actually reduce the total dispersion by 

Qs ' delaying the formation of a concentration shock. However, if the electro-osmotic flow is 

sufficiently high, the total dispersion is increased because of the Taylor- Aris contribution. 

cn 

o 

psj ■ 1. Introduction 

Macromolecules such as DNA and proteins carry an electrical charge in aqueous so- 
lution because of ionic dissociation of molecular groups on their surface. The amount 
'k^ i and sign of the charge depends on various factors such as the pH of the solution. If an 

'V^ I external electric field is applied, these macro-ions migrate along the field with a velocity 

C^ i that depends on its charge, size and shape as well as on the ionic composition of the back- 

ground electrolyte. Electrophoresis is a technique of analytical chemistry for separating 
a mixture of chemical compounds by exploiting their differing migration velocities in an 
applied electric field. It is a widely used laboratory technique in areas such as molecular 
biology, forensics and medicine. 

The most common format is "gel electrophoresis," where the electrolyte is suspended 
in a porous gel. The technique of capillary electrophoresis (CE) has developed rapidly 
in recent years, partly because of the possibility of integrating CE channels on to micro- 
fluidic chips. In CE, the sample and suspending electrolyte are contained in a micro-fluidic 
channel with width in the range of tens of microns. CE can proceed in several modes, the 
simplest of which is "zone electrophoresis" . Here a sample zone or band is injected into a 
micro-channel, which then moves by electrophoresis towards the detector at the opposite 
end under an applied electric field. Different components arrive at different times and 
their arrival is marked by a peak in some solute property (usually the UV absorbance) 
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at the detector window ( Landerall996tlCamillerilll998l ). The wahs of the capillary, which 
are most commonly made of fused silica, has an electrostatic charge that is characterized 
by its zeta potential. Th erefore, the appU ed electric field results in an electro-osmotic 
flow along the capillary (jProbsteinI 119941 ) . This flow is normally advantageous in CE 
because it sweeps both positive and negative ions past a single detector near the outlet 
and also reduces the transit time from inlet to detector. However, it can, under certain 
circumstances adversely affect the separation by causing "anomalous dispersion" of the 
sample peak (jGhosal 2006) . 

The physical process of separation in CE is the result of the mutually opposed and 
competing processes of differential migration of ions and diffusive spreading along the 
axis of separation. The resolution depends on the strength of the applied electric field. 
In practice, the applied voltage is often very large, in the range of kilo volts. The highest 
voltage that can be applied is limited by the ability to dissipate the considerable Joule 
heat that is produced in the buffer by the electrolytic current. This is the reason why 
the channel width may not exceed some tens of microns. It is clearly advantageous to 
have an electrolytic buffer of low conductivity in order to minimize heating. It is also 
advantageous to have a relatively high concentration of sample ions, since one of the 
limitations of CE is the high demand placed on the sensitivity of the detector, which 
must be sensitive to light attenuation over an optical path length of only a few tens 
of microns. However, the ratio of sample to background ion concentration cannot be 
increased indefinitely; distortions due to "electromigration dispersion" or the "sample 
overloading effect" limits the highest ratio of sample to background ion concentration 
that may be used in CE. 

The underlying ph ysical mechanism of electromigration dispersion was elucidated by 
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Mikkers et al\ (|19790 and may be roughly explained as follows: when the concentration 
of sample ions is comparable to that of the carrier electrolyte, the local electrical conduc- 
tivity is altered in the vicinity of the sample peak. On the other hand, since the current 
through the capillary must be the same at all points along the axis, the electric field must 
change locally. This follows since the product of the conductivity and electric field is the 
current (according to Ohm's law, neglecting for the time being the diffusion current due 
to ionic concentration gradients). This varying electric field alters the effective migration 
speed of the sample ions, which, in turn, alters its concentration distribution. Thus, we 
have a nonlinear transport problem that must be solved in a self consistent manner. 
Electromigration dispersion causes highly asymmetric concentration profiles, high rates 
of dispersion and shock like structures that are reminiscent of nonlinear waves seen in 
many other physic al systems. These effe cts have been widely reported in the literature 
on electrophoresis (|Bouskova et al\\2004] . 

Simple o ne dimensional mathematical models o f electromigration disp ersion have been 
studied bv iMikkers et al\ (|l979l ): iMikkerd (|l999l ): iBabskn et ail (|l989( ) by invoking the 
assumption of vanishing diffusivity which results in a single nonlinear hyperbolic equa- 
tion for the concentration of sample ions. Solutions describe the observed steepening of 
an intit i ally sm ooth profile leading to subse quent shock formation. Two recent reviews 
by ICaa ([20091) and iThormann et ali (J2009|) provide more exte nsive references t o this 
early work. The restriction of zero diffusivity was removed bv IChosal fc ChenI ( 20101 1 
(henceforth referred to as GC) who considered a three ion system - the sample ion, co- 
ion and counter-ion - where the diffusivities of the three ionic species were equal but 
not necessarily zero. The sample concentration was shown to obey a one dimensional 
nonlinear advection diffusion equation which reduced to Burgers' equation if the sample 
concentration was not too high relative to that of the background ions. Since the ini- 
tial value problem for Burgers' equation may be exactly solved, useful insight into the 
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Figure 1. Schematic diagram (not to scale) illustrating the mathematical problem of electro- 
migration dispersion of a sample in CE. In addition to the sample ions shown, the capillary also 
contains a carrier electrolyte consisting of co- and counter- ions. 



nature of electromigration dispersion could be gained in the limits of small as well as 
large Peclet numbers. A generalization of this model to the situation where the buffer is 
a weak electrolyte has recently been presented (jChen &: Ghosalll2011l ) . It was found that 
the evolution of the peak shape may still be described by the Burgers' equation with 
a shght re-interpretation of a certain model parameter. If this parameter is taken as a 
fitting parameter, excellent quantitative agreement is obtained with measured values of 
the peak v ariance in experime nts. A similar nonlinear advection diffusion equation was 
derived bv lZangle et all ( 20091 ) in order to describe "desalination shocks" near constric- 
tions in micro channels. However, in their problem the physical origin of the nonlinearity 
is related to surface conduction effects. 

In this paper, we extend the earlier work of GC by taking into account the fact that the 
channel wall generally has a non-zero zeta potential which creates a bulk fluid flow in the 
capillary due to electro-osmosis. The electro-osmotic flow affects the transport process 
in the following way: the axial variation of the electric field created by the conductivity 
changes results in a variable electro-osmotic slip velocity on the channel walls. This in 
turn results in an induced pressure gradient and the concomitant appearanc e of a shear 



in the ve locity profile across the channel due to a well known mechanism (jHerr et 



2000t iGhosal ^ 2002 a. 6 . c) . The shear then enhances axial dispersion due to the Taylor- 
Aris effect ( Tavloiill953t lArisI 119561) . Our principal result is equation p.23p . which is 
a one dimensional nonlinear advection diffusion equation for the sample concentration 
averaged over the channel cross-section. 

The rest of the paper is organized as follows. In § [2] the complete mathematical state- 
ment of the flow and electro-diffusion problem is presented. In §|3]a systematic reduction 
of the full equations leading to equation p.23p is achieved by introducing a number of 
physical approximations. The consequences of this equation are discussed in § |4l where 
it is shown that there is a peak in the efficiency of separation at an intermediate value of 
the electro-osmotic flow strength. The broad implications of the analysis to separation 
efficiency in CE are discussed in § [5l 



2. Mathematical Formulation 

We consider (see Figure 1) a uniform, long capillary of arbitrary cross-sectional shape 
connected to infinite reservoirs at either end. The length of the capillary is much greater 
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than the width; the capillary may, for most purposes, be assumed infinitely long. The 
capillary contains a solute and N strong (i.e. completely dissociated) electrolytes of con- 
centration Ci{x, J/, z, t) where i = Ito N. These concentrati ons then obey th e conservation 
equations based on the Nernst-Planck model for ion flux ( Probsteinlll994l ) 

uc 

— ^ + V • (uc, - z,eiyiCiV<j)e - AVc,) = 0, (2.1) 

ot 

where Zi, Vi and Di are the valence, mobility and diffusivity of species i, the electronic 
charge is e, i^e is the electric potential and u(a;,y, z,t) is the fluid flow velocity. The 
potential 0e obeys the Poisson equation (in CGS units) 

N 

1=1 

where e is the permittivity of the solvent and pe is the density of free charges in the 
solvent. Since the Reynolds number is small for flow in micro-capillaries, the fluid flow is 
described by Stokes equation 

-Vp + T]V^U-p,V^e=0, (2.3) 

together with the incompressibility constraint 

V-u = 0, (2.4) 

where p is the pressure and rj is the dynamic viscosity of the solvent. 

These equations are associated with a number of boundary conditions that describe 
our situation. The concentrations Ci reduce to their equilibrium distributions in a uniform 
capillary very far away from the sample zone 

Ci(a;,y,z) ~cf(y,2;) if |a::| ^ oo. (2.5) 

Clearly c1'^(jj, z) — when i corresponds to the sample and is given by solutions of the two 
dimensional Poisson-Boltzmann model for the other species (background electrolytes). 
The potential 4>e obeys Dirichlet boundary conditions 

{0^ if (x, y, z) e Snj 
V a {x,y,z)e S+ (2.6) 

ii{x,y,z)eS- 

where (j>'^ is the value of the potential at the wall, V is the applied voltage, and, S^, 
5-1- and S- are the parts of the bounding surface of the capillary that correspond to 
the walls, the inlet and the outlet respectively (Figure 1). The inlet and the outlet 
sections are separated by a distance L and we are interested in the limit where V and 
L approach infinity but V/L = Eq remains finite. The hydrodynamic flow satisfies the 
no-slip conditions at the wall 

u{x,y,z,t)^0 ii {x,y,z)eSiu, (2.7) 

and, 

u{x,y,z,t)^u''''{y,z) if |2;| ^ oo, (2.8) 

where u'^'(y, z) is the electro-osmotic flow in the uniform capillary when no sample is 
present, and, when L is inflnite. Finally, the ion fluxes are zero at the capillary walls 

(z^ei^jCiV^e + AVci) • n = ii {x,y,z)e S^ (2.9) 
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where n is the unit normal at the wall directed into the fluid. 

These equations and boundary conditions form a well posed system of equations that 
may be integrated numerically for given initial conditions. Such a description would, how- 
ever, be unnecessarily complex. Our objective is to introduce a series of approximations 
that enable a reduced description in terms of a single one dimensional evolution equation 
for the cross-sectionally averaged concentrations Ci(x, i), since this is the quantity that is 
essentially measured by the detector. Henceforth, a bar above a field variable will always 
denote its average over the cross-section of the capillary. 



3. One dimensional homogenized equations 

We now introduce a number of approximations that simplify the description and lead 
up to the one dimensional homogenized equation for the sample ion concentration that 
we seek. The developments closely follow GC except for the important difference that 
the possibility of a hydrodynamic flow field u is admitted. 

3.1. Thin Debye layers 

The fixed charges on the wall are shielded by a cloud of counter-ions over a length scale 
of the order of the Deb ye length (A^) which is related to the equilibrium ionic concentra- 
tions ( Probsteinlll994r ). The Debye length is typically of the order of nanometers; much 



smaller than the characteristic width of the micro-channel which is in the range of tens 



of microns. Under such conditions, the left hand side of (|2.2[) may be set to zero (jPlanck 
Il890 1. and, as a consequence, the last term in (12. 3p . representing the density of electric 
forces also vanish. Electrical forces coul d also arise outside of Debye layers in res ponse 



to variations in electrical conductivity ( Chen et aLll2005t lOddv fc Santiagol 120051 ). The 



contribution from such forces may however be neglected in the current context (see 
Supplementary Material online). Thus, Poisson's equation is replaced by the constraint 
of local electro-neutrality. The loss of the equation for (pe is however not catastrophic, 
because ^e can be determined from the equation of current conservation 



V-Je = (3.1) 



where the current density 



N \ N 



Y^ zle^ViC^ V(/)e - ^ z.eA Vc, (3.2) 



is the sum of the Ohmic conduction (first term) and a contribution due to differen- 
tial diffusion (second term) . Equation p.ip follows from the constraint of local electro- 
neutrality on summing (|2.1[) after multiplying by z^e. The loss of the electrical forcing 
term in (j2.3p is mitigated by replacing the no-sli p boundary con dition (j2.7p with the 



Helmholtz-Smoluchowski slip boundary condition (|ProbsteirJll994l ) 



u(x,j/,z) = — ]i{x,y,z)£Sw (3.3) 

Airrj 

that takes into account the presence of a boundary layer of non-vanishing charge next 
to the capillary wall. Here C, the so called "zeta-potential" is the difference in the values 
of (pe between a point at the outer edge of the Debye layer and the corresponding point 
on the wall. It is a material property that will be assumed constant. The boundary 
conditions (12.61) are replaced by 

Je • n = (3.4) 
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which follows from ()2.9|) and is the condition of zero current flow into capillary walls. At 
large distances from the sample zone we have a uniform field (-Bo) and a uniform flow 

(Ueo), thus, 

(j)e{x,y, z,t) ^ —Eqx if |a;| -^ oo (3-5) 

and 

u(x, y, z, t) ^ u'^'^{y, z) — x = Ueo x if |a:| — >■ oo (3-6) 

where x is the unit vector in the axial direction. 

3.2. Three ion model 



We introduce the Kohlrausch ( Kohlrauscblll897l ) function 



^ 1 
K{x,y,z,t) :^S2—Ci{x,y,z,t). (3.7) 

2—1 

Then (|2.ip together with the constraint of local electro-neutrality yields the following 
equation for K: 

For a two ion system, equati on (13.81) together with the condition of local electro-neutrality 



leads to the Ohmic model ( Melcher fc TavloiJll969t IChen et al\\2003t \ where the local 



electrical conductivity evolves as a passive scalar with an effective diffusivity D^ = 
2D^D_/{D^ + D^); D^ and I?_ denotes the diffusivities of the cation and anion. In our 
problem, we have more than two ionic species and the Ohmic model is not applicable. 
However, if we assume that the ions all have the same mobilities, Vi = v, and therefore, 
on account of the Einstein relation Di/vi = ksT, the same diffusivities, Di = D, then 
the function K , rather than the conductivity, is found to evolve as a passive scalar: 

V + V • (uK) = DV^K. (3.9) 

ot 

We restrict ourselves to a minimal model consisting of a system of just three ions {N = 3). 
We will drop the index i and instead use the suffix p and n to denote the positive ion 
and negative ion respectively, in the background electrolyte. The absence of a suffix will 
indicate the sample ion. For example, Cp, c„ and c are the concentrations of positive 
ions, negative ions and sample ions respectively. Then from the local electro-neutrality 
constraint and the definition of K we have 

c„ _ ct') = I^^c + ^^^^ SK, (3.10) 

c^ + ^ci^) = ^l^Ij^c~^^^^SK, (3.11) 

Zp Zp Zji Zp Zji 

where Cn is the concentration of negative ions in the background electrolyte far from 
the sample zone. If K^ denotes the value of K in the far field, then the perturbation 
6K = K — Kao is advected by the flow and spreads diffusively from the injection region. 
Therefore (GC), after an initial transient that is small compared to the total analysis 
time, the sample peak, which moves by electrophoresis in addition to being advected 
by the flow, migrates into a region of space where SK is effectively zero. Thus, in the 
vicinity of the sample peak, we may assume that SK = 0, so that, the background 



Electromigration Dispersion 7 

ion concentrations Cp and c„ may be expressed as linear functions of the sample ion 
concentration c. 

The diffusion current represented by the second term in p.2|) vanishes when Di = D 
on account of local electro-neutrality, so that it reduces to Ohm's law for a homogeneous 
electrolyte: Je = — CeV^e, where 



N 



2:^6 i^iCi = (Too (1 - a^) (3-12) 

4=1 



is the local electrolyte conductivity. The expression on the right is obtained on elimi- 
nating Cp and c„ using p. 101) and p. Ill) with 6K = 0. Here tJoo is the bulk electrolyte 
conductivity, (p — c/cn is the sample concentration relative to the background and a 
is the parameter 

(z — Zr,)(z — Zri) 

a=^ ^iZ^ — ^ (3;L3) 

^n \^p ^n) 

introduced in GC to characterize the nature of the nonlinearity. 

3.3. The lubrication limit 

The subsequent development is based on the premise of "slow axial variations" on account 
of the inlet to detector separation, L, being much larger than the characteristic channel 
width, Wo- In practice, L is of the order of tens of cm whereas wq ^ 10 — 100 fj,m. Thus, 
L/wq ~ 10^ — 10"*. Therefore, as the sample moves down the micro-channel, the sample 
concentration and all other quantities controlled by this concentration vary on an axial 
length scale L^ 3> wq. To see this, suppose that the initial sample concentration was a 
delta function in the axial direction. Then it would spread over a distance of the order 
of the channel width in time t^ '^ Wq/D. The total analysis time is Tq ~ L/vq, where vq 
is a characteristic migration velocity. Thus, Ta/rd ~ (L/wo)Pe~^ where Pe — vqWo/D is 
the Peclet number. Since typically Pe ~ 10 — 100, Ta/rd ;» 1 so that concentrations are 
homogenized across the channel and the fluid flow and transport problems both become 
quasi one dimensional a short time {'^ Td) after sample injection. 

The equation of current conservation, p.ip . may then be integrated using the boundary 
condition p.4|) to give 

E = x£;o/(l-a0) + ---, (3.14) 

where E = ~'^4'e is the electric field and the neglected terms are asymptotically small in 
the lubrication limit. Thus, the field is predominantly in the axial direction and depends 
on the local sample concentration. The axial variability of the electric field affects the 
evolution of the sample concentration in two ways. First, it results in a variable elec- 
trophoretic migration velocity wo/(l — «</') for the sample ions, where vo = zcuEq is the 
migration velocity of an isolated sample ion. Second, it results in a variable slip velocity 
for the electro-osmotic flow through the boundary condition p.3p . 

The hydrodynamic fio w field due to su ch a variable slip velocity may be written down 
using lubrication theory (jGhosall 12002 cl ). The velocity is predominantly in the axial di- 



rection, u = MX + • • • and the axial component u may be written as the sum of a mean u 
and a fluctuation about the mean, Au, due to the induced pressure gradient. The mean 
flow -0 is a constant given by 

1 [^ 

u— lim — / Ug{x,t) dx = u^o (3.15) 

L^oo L Jq 

where Us{x,t) is the slip velocity at the wall. In the presence of concentration gradi- 
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ents, the slip velocity could also have a diffusiophoretic component ( Prieve et aLlll984i 
Rica fc Bazantll2010l ). This is however expected to be small in the present context (see 
Supplementary Material online). 



Us{x,t) = Ueo/(l 



and, 



Au — U — U = {Ug — Ue 



■'|i-f1- 



1 — a- 






+ 



(3.16) 



(3.17) 



to leading order in the lubrication approximation. Here Up is a function that depends 
solely on the cross-sectional shape and represents the flow profile due to a unit pressure 
gradient. It is defined by the solution of the equation 



9^M<o d'^u„ 



dy"^ 



dz 



= -1 



(3.18) 



with the boundary condition 



Up(y, z)=0 \i{x,y,z)e Sw (3.19) 

The function Up admits analytical representation for several cross-sectional shapes. 

3.4. Taylor-Aris limit and macrotransport 
The time evolution of the concentration field of an advected scala r in the limit Ta/Tij ^ 1 



was first described by Taylor and Aris (|TavloiJll953l : lAriall956l) and later applied to a 
wide variety of problems involving s hear induced dispersion ( Brenner &: Edwardall993l) 
including dispersion problems in CE ( Ghosalll2006l ). In this limit, lateral inhomogeneities 
in the scalar concentration are small c{x, y, z, t) = c(x, t)+- ■ •, so that the cross-sectionally 
averaged concentration c is advected by the mean fiow u and undergoes axial diffusion 
with an effective diffusivity ( Datta fc Ghosalll2008l ) 



D,=D 



Gu 

IT 



where G satisfies 



and the conditions 



d^G d^G 



dy^ 



dz 



= Au = 



(3.20) 



(3.21) 



VG-n = ii{x,y,z)eS^,. (3.22) 

To remove the indeterminacy of G up to a constant, we further impose the condition G = 
0. Furthermore, the sample ions are advected with a concentration dependent velocity 
^o/(l ~ acf)). Thus, the equation satisfied by (f) is 



dt 



d_ 
dx 



«o 



1 — a0 



d_ 
dx 




(3.23) 



where fc is a numerical constant depending solely on the cross-sectional shape and wq is a 
characteristic channel width. The expression for the effective diffusivity D^ appearing on 
the right hand side of (|3.23p follows on evaluating p.20p using (|3.17p . In particular, for a 
planar channel of half-width wq, a simple calculation shows that Up = (wg — y^)/2, which 
leads to fc = 2/105. Equation p.23p is the generalization of the evolution equation derived 
in GC to the situation where the channel has a non-zero electro-osmotic slip velocity. 
An alternative derivation of p.23|) using the method of multiple scales is sketched in the 
Appendix. 
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Figure 2. Time evolution of the normalized rate of increase of variance for three different 
values of the dimensionless electro-osmotic flow strength u« = u^o/vo- Horizontal dotted line is 
the effective diffusivity predicted by equation (35) in GC. 



If \a\(f) 4; 1, the nonlinear terms in (|3.23p may be expanded in Taylor series: (1 — 
a4>)~^ = 1 + a^ + a^cjr + • • •, so that, in place of Burgers' equation arrived at in GC we 
get the following equation 



'dt 



+ {vo + Ueo + 2avo<j) + 3a'^vo(f>^) 



d(t)_ 
dx 



d_ 
dx 



D 



D ' 



dx 



(3.24) 



which shows an amplitude dependent contribution to the diffusivity in addition to a 
correction to the term 2avo4>{d(f)/dx) representing nonlinear wave steepening. Equation 
(|3.23p . or, its weakly nonlinear version p.24p . is our principal result. 

Equation (|3.23p has a singularity at (j) — 1/a when a is positive. However, (|3.23l) ceases 
to be valid even before this singularity is reached (GC), because, the requirement that Cp 
and c„ must be non-negative, imposes the constraint (j) < (p^ where 4>c = {zp — Zn) / {zp — z) 
when z < and (f>c = — [z„(zp — z„)/[zp(z — z„)] when z > 0. If a > 0, then it may be 
shown that (pc < 1/a. The physical reason for the breakdown of p. 231) when cj) > (j>c is 
the following: under those conditions the conductivity in the sample zone is so high as 
to reduce the electric field to very small values. Then the electrophoretic motion of the 
sample ions become very small and there is little relative motion between the sample 
peak and SK. Thus, the assumption that 6K may be set to zero in the sample zone 
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can no longer be employed and p.23p is no longer valid. However, after sufficient time 
has passed and axial diffusion has caused the peak value of (f> to fall below 0c, the time 
evolution once again proceeds in accordance with (J3.23I) . Numerical simulations of the full 
electrohydrodynamic equations described in i; 13. II confirm this behavior (Chen & Ghosal, 
unpublished) . For the purpose of numerical integration, (I3.24p is a little more convenient 
than p.23p since it does not exhibit the singularity at (f> — 1/a. Nevertheless, even though 
solutions to p. 241) may be formally calculated even for (j) > (pc, the solution in this 
regime is devoid of physical significance, and indeed, will yield negative concentrations 
of background electrolytes in parts of the domain if (I3.10p and (13.111) with 6K = are 
employed to calculate the concentrations of the background ions. 



4. Dispersion 

Equation (13.231) or (j3.24p provide a compact description of electromigration dispersion 
that is helpful for gaining a qualitative understanding of the underlying mechanisms. 
Furthermore, it is much more amenable to numerical integration than the full three 
dimensional coupled problem involving fluid flow and transport. By way of example, in 
this section, we use numerical integration to illustrate a point that at first sight may 
appear counter-intuitive, but, is clarified by an analysis of (|3.23|) . 

Since the presence of a zeta-potential results in cross-channel variations in the flow 
velocity due to induced pressure gradients, one would expect the efficiency of separation to 
be adversely affected. However, in certain ranges of parameters the reverse may actually 
be true. To understand this, one needs to examine the roles of the different terms in 
p.24p . The nonlinear term on the left hand side causes wave steepening leading to the 
formation of shock like structures. This is the dominant mechanism that contributes to 
electromigration dispersion (GC). Taylor dispersion can actually mitigate this tendency 
by increasing the effective axial diffusivity that serves to diffuse the electrokinetic shock. 
However, on the other hand, if the Taylor dispersion is too large, its contribution to the 
axial dispersion dominates with a consequent loss of separation efficiency. Thus, there 
is an intermediate value of the wall zeta potential that corresponds to the lowest peak 
dispersion. 

Equation (j3.23p was integrated numerically using a finite volume method that allows 
for adaptiv e grid refinement and varia ble time steps. We used the "odelSs" solver in 
MATLAB ( Shampine fc Reicheltlll997l) (see Supplementary Materials online for further 



details). The geometry chosen was that of a planar channel of width 2wo- A "moving 
window" that is advected with the mean migration velocity vq + u^o was used to optimize 
the number of grid points needed. The size of the window was chosen large enough that 
(j) can be set to zero at the domain boundaries with negligible loss in accuracy. The initial 
concentration profile 0(a;, 0) was chosen as a Gaussian of standard deviation wq centered 
on a; = and with different peak strengths. The "sample loading" is characterized (GC) 
by a Peclet number P — vqT/D based on the length scale 

r= / (t)dx, (4.1) 



which is an integral of motion. A second Peclet number that characterizes the diffusion, 
Pe = vqWq/D is held fixed at the value 200. The parameter a was set to 0.5. These 
conditions are fairly typical for laboratory experiments. 

Figure [2] shows the evolution of the quantity Dt = {2D)~^da^ /dt, where cr^ is the 
variance of (/)(a:;, t), as a function of the dimensionless time vqI/wq. At late times, when (/> 
is sufficiently small, p.23p may be replaced by p.24p . At even larger times, the quadratic 
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Figure 3. The number of theoretical plates N — I? ja^ as a function of dimensionless elec- 
tro-osmotic velocity u, = Ueo/vQ for three different detector to injection point distances (L) 
and a range of sample loading characterized by the Peclet number P = vqT / D. Here A^o is the 
"ideal" value of A'^ in the absence of electromigration dispersion. 



terms in become vanishingly small, so that, the equation describing the evolution of 
(j) essentially reduces to Burgers' equation as discussed in GC with the minor difference 
that the constant part of the advection velocity is uq -|- Ueo and not vq. Therefore, one 
would expect that Dt, which is like an "instantaneous" difFusivity, should asymptote to 
the value given by equation (35) in GC. Figure [2] shows that indeed, this is the case 
for all values of the flow strength, m, = Ueo/wg. However, the evolution of Dt to the 
common asymptotic value follows different trajectories depending on the strength of the 
electro-osmotic flow. For larger values of w* , Dt is initially relatively large because of the 
contribution of the quadratic (and higher order) terms in (j) to the effective diffusivity in 
(|3.23p . This is because of Taylor dispersion. At intermediate times, Dt is lower than the 
asymptotic value predicted by equation (35) in GC, a consequence of the fact that a larger 
effective axial difFusivity delays the formation of shock like structures that result from 
nonlinear wave steepening. The total variance is the initial variance (wq) plus the area 
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Figure 4. Concentration profiles <j>{x, t) at a fixed instant of time {vot/wo = 200) for several 
values of the electro-osmotic flow strength, u* = Ueo/vo. The increased effective axial diffusivity 
due to Taylor dispersion "softens" the electromigration shock that tends to form at the leading 
edge. Here P = 50 and a; = Xc is the location of the centroid of the peak. 



under the curve Dt from t = to some time t — tf when the peak arrives at a detector 
located at x — L. This is shown in Figure 2 and 3 in the Supplementary Material. 

The efficiency of separation in CE is often characterized by the "number of theoretical 
plates" N — L"^ /a"^ which is a dimensionless measure of the peak width. In the top panel 
of Figure 131 we plot iV as a function of w, for a number of downstream detector positions 
and for sample loading ranging from P = 0.5 (weak) to P = 50 (strong). It is clear that 
increasing P reduces N due to electromigration dispersion. However, the dependance on 
u* is non-monotonic with N increasing at first with u* but then decreasing or leveling 
off to a plateau depending on the location of the detector. Increasing the electro-osmotic 
flow Ueo generally increases N due to an obvious and quite trivial reason. If the detector 
position is fixed, the peak reaches it faster for higher electro-osmotic flow speeds and 
the accumulated variance is less simply because the peak has evolved for a shorter time 
period. In fact, in the absence of the nonlinearities induced by electromigration dispersion, 
0-2 EE ctq = Wq + '^Dtf ^wl + 2DL/{ueo + Vo) ■ In this ideal limit N = N^ ^ L^/^^l- In 
order to eliminate this obvious effect that electro-osmotic flow has on peak dispersion, 
we plot N/Nq as a function of u* in the lower panel of Figure |31 which clearly shows, 
that the degradation of separation efficiency by electromigration dispersion is minimized 
at an intermediate value of the electro-osmotic flow w,. The dependence of A^ on L at 
different values of P is shown in Figure 4 of the Supplementary Material included in the 
online version. 

The mechanism of this reduction in total variance when u* is sufficiently small is evi- 
dent from Figure |4] depicting peak shapes at a fixed time Vot/wo = 200 for a number of 
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Figure 5. Concentration (f){L,t) as a function of time of arrival (f) at a fixed detector location 
{x = L) for several values of electro-osmotic flow strength, u* — Ueo/vo- It is seen that though 
electro-osmotic flow improves the resolution, the degree of improvement "saturates" as the flow 
rate is increased, as one might expect from Figure [S] Here P = 50. 



values of u*. It is seen that the increase in the effective axial diffusivity due to Taylor 
dispersion "softens" the electrokinetic shock that results from nonlinear wave steepen- 
ing due to electromigration. Thus, electro-osmotic flow has three competing effects (a) 
reduction in the time available for diffusion (b) reduction in the nonlinear wave steepen- 
ing on account of enhanced effective diffusivity (c) increase in axial diffusion due to the 
Taylor- Aris effect. The total dispersion is determined by the relative contribution of each 
of (a), (b) and (c). At low values of m*, (a) and (b) dominate, but at higher values, (c) is 
more important. In laboratory systems, u* could have either sign and a magnitude that 
could vary between zero (e.g. in a coated capillary where electro-osmosis is suppressed) 
to some number of order unity, since the electro-osmotic velocity is similar in magnitude 
to the electrophoretic velocity. Figure [5] provides an alternate viewpoint that has a closer 
correspondence with experiments: the signal intensity (t>{L, t) is plotted as a function of 
the time of arrival at a fixed detector position at distance L from the injection point. It is 
seen that introduction of an electro-osmotic flow at flrst results in a significant reduction 
in peak width, but with diminishing returns as the fiow strength is increased. A similar 
plot is shown as Figure 6 in the Supplementary Material except there (^(x, t) is plotted 
as a function of x at fixed times. 



14 S. Ghosal and Z. Chen 

5. Conclusions 

In conclusion, we would like to make a few remarks about the relevance of the analysis 
presented here to laboratory practice. The model studied here is clearly oversimplified 
and (|3.23p cannot be expected to predict in detail the peak shapes in a real laboratory 
experiment. In particular, real electrophoresis buffers contain many more than three ions 
including one or more weak acids or bases that are added to maintain a stable pH. Nev- 
ertheless, it is important to understand the qualitative effects predicted by the simplified 
model considered here, since the actual dispersion in a real laboratory experiment is likely 
to be due to a multiplicity of causes, of which, the mechanism discussed here would be 
one. 

In this paper, we have described the effects of electromigration dispersion by means of 
a small set of dimensionless parameters: the Peclet numbers Pe = vwq/ D, P — vT/ D, a, 
and, the dimensionless electro-osmotic fiow speed u* = u^o/vq. A rough idea of typical 
values of these parameters in laboratory experiments is helpful. If we take vg '^ 1 mm/s, 
Wq ^ 10/im, D ^ 10~^ cm^/s as fairly typical, we have Pe ^ 10. For large molecules, 
D could be an order of magnitude smaller, so that Pe could be in the hundreds. If the 
sample concentration is of the same order as the background electrolyte concentration 
and the peak width is of the order of the channel width, then F ~ wq. However, such 
sample concentrations would be considered extremely high in CE and would show very 
strong peak distortion. We may therefore take this as an upper limit for F. Thus, P 
could range from essentially zero (the linear regime) to P = vT / D ~ vwq/D — Pe. Since 
electrophoresis and electro-osmosis are generally of comparable magnitudes, u^o '^ v. 
Thus, w* '-^ 1, though u* could be close to zero if the channel walls have a low zeta 
potential (as in a coated capillary). The magnitude of the parameter a is of order unity if 
the sample valence and the valence of the background electrolytes are similar. However, 
for macro-ions \z\ ~ 10 or larger. For such large value of z, a ~ z^ could be in the 
hundreds, and, nonlinear effects of the kind considered in this paper would be very 
pronounced. 

Electro-osmotic flow is generally desirable in the context of CE as it shortens the transit 
time from inlet to detector and as a consequence also reduces the peak variance. The effect 
of electro-osmotic flow on separation efficiency is however complicated in the presence 
of electromigration dispersion since a number of competing effects are simultaneously 
present. This paper generalizes the analysis presented in an earlier paper (GC) where 
the effect of electromigration dispersion in the absence of any bulk flow was considered. 
However, capillary electrophoresis channels are usually made of fused silica which has a 
large negative zeta potential and sustains strong electro-osmotic flow. Thus, the theory 
presented here broadens the scope of the earlier paper making it more relevant to practical 
systems. 

In summary, a homogenized equation for the cross-sectionally averaged concentration 
of sample ions was derived by exploiting the approximations outlined in § [3l which are: 
(a) the Debye length is much smaller than the typical channel width (b) the electrolyte 
contains only three ionic species of equal diffusivity which are also strong electrolytes 
(c) the time between injection and detection is much longer than characteristic trans- 
port times across the capillary. The outcome of the analysis is the replacement of the 
molecular diffusion coefficient in the one dimensional transport equation derived in GC 
by an effective axial diffusivity which is a nonlinear function of the concentration. If 
the concentration is not too high, this effective diffusivity is a quadratic function of the 
concentration. 

The reduced one-dimensional equation p.23|) provides a compact description of elec- 
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tromigration dispersion that is much more amenable to numerical integration than the 
full coupled three dimensional problem involving fluid flow and transport in the capil- 
lary. In fact, since both the length to width aspect ratio of the capillary as well as the 
Peclet number (Pe) are usually large in experiments, such three dimensional simulations 
can be computationally intensive. The loss of accuracy in replacing the full equations by 
the one dimensional model is minimal because the requirements for its validity are well 
satisfied and the asymptotic convergence of the thin Debye layer and long length scale 



approximations have been shown to be rapid ((GhosalluOOGt iDatta fc: Ghosal 2009. 2008 ) 



3Dye 

Si 



The effect of a wall zeta potential may be explained qualitatively in the following 
way: alteration of the ionic composition of the solute in the sample zone results in a 
change of the electrical conductivity and therefore in the axial electric field. Thus, the 
electro-osmotic slip velocity now has an axial variation which induces pressure gradients 
along the channel. This in turn induces shear in the electro-osmotic velocity profile. The 
consequent shear induced (Taylor- Aris) dispersion results in a concentration dependent 
axial diffusivity. Numerical integration of the homogenized equation p.23p provided some 
qualitative insights into the effect of this nonlinear diffusion on separation efficiency. It 
is seen that the total variance for a fixed detector position actually decreases with the 
strength of the electro-osmotic fiow as long as the fiow is not too strong. This is due to 
the fact that an increase in the effective diffusivity for axial transport counteracts the 
wave steepening resulting from electromigration dispersion, and, furthermore, the transit 
time from injection to detection is reduced. However, if the ffow is much stronger, the 
dispersion caused by the enhanced axial diffusivity itself is dominant, and this degrades 
any gain in resolution due to the aforemetioned causes. Thus, if all other parameters 
are invariant, the peak variance is a non-monotonic function of the electro-osmotic fiow 
strength and is optimal at a certain intermediate value of the ffow velocity. 

This work was supported by the NIH under grant R01EB007596 
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Appendix A. Derivation of (13.23^ by the method of multiple scales 

It is advantageous to work with dimensionless variables. This is achieved by using wq as 
the unit of length and wq/vq as the unit of time. The governing equations are therefore 
those presented in §[2] and §[3] but with wq and vq set to one. These equations are 

[(1 - acj))Xl + [(I - a(b)Y]y + [(1 - a0)Z], =0 (A I) 

(l>t + [{U + X)(P],^ + [{V + Y)(P]y + [{W + Z)(f>], = Pe-\(f)^^ + (l)yy + (j) , ^] (A 2) 

-PX + UxX + Uyy + U^2 = (A 3) 

-Py + Vxx + Vyy + Wz2 = (A 4) 

-PZ + Wxx + Wyy -f W^^ = (A 5) 
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where {X,Y,Z) is the electric field vector (normalized by £^o), u — {u,v,w), (j) is the 
normalized sample concentration, Pe = vqwq/D (Peclet number) and suffixes denote 
partial derivatives. The flow is solenoidal: 

Ux + Vy + Wz — Q (A 6) 

and the electric field is irrotational: 

Xy = Yx Xz = Zx Yz = Zy (A 7) 

The boundary conditions at the wall are 

u^u^X v = w = (A 8) 

and 

mcpy + n(j)z 
mY + nZ 

where (0, m, n) is the unit normal on the wall and u 
and downstream of the sample peak 

X-1 u-u* w-0 w-0 (/)-0. (All) 

This description presumes infinitely thin Debye layers (§ 13. ip and constancy of the 






(A 9) 





(A 10) 


i W* : 


- Ueo/vQ. At points far upstream 



Kohlrausch function in a three ion system (ii l3.2[) . 

We now assume that all axial variations happen on a slow spatial scale for which we 
introduce the slow variable £, = ex and the corresponding slow times t = et and T — e^t, 
where e is a small parameter. The physical origin of the two times is due to the existence 
of the time-scale wq/vq which is the advective time over a channel diameter and the 
diffusion time-scale Wq/D for diffusion across the channel diameter. Since the flow and 
fields are predominantly axial, we introduce 

Y = Y/e Z = Z/e ii ^ v/e w^w/e p^ep (A 12) 

where the pressure scaling is designed to retain the pressure gradient term at leading 
order. Next we re- write (jA ip - (jA5P in terms of the scaled variables using the transforma- 
tions dt = edr + c^Ot and dx = e9j and expand all dependent variables / in asymptotic 
series / = /^''•' + e/^^-* + e^/*-^' + • • • Equating like powers of e on cither side of the 
governing equations we get a series of successive problems to solve. 

A.l. Zeroth Order 
From (JA1|1 and (|A7jl 

[(1 - a0("))x(°)]e + [(1 - a0(°))y(")], + [(1 - a0(°))z(o)], = 0, (A 13) 

Xf ) = ^(0) = 0. (A 14) 

Thus, X'") is independent of y and z. Integrating (jA 13[) over the channel cross-section 
and using the boundary condition (|A 10|1 we get 

[(l-a0(°))x(o)]^ = O, (A 15) 

and thus, on account of the far field conditions (jA lip 

X(") = l/(l-a^(o)). (A 16) 

Equation (JA2| gives at lowest order 

</'^°J+</'£'=0, (A 17) 
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which, with the Neuman boundary condition (jA9|) impUes that (j)^ ' is independent of y 
and z, thus, 

0(")=0(o)(e,r,T). (A 18) 

Finally, the flow equations (|A3P - (jA5|) give 



6r+<^+4?=0, (A 19) 



Pi 



-^r^-^^ -0- (A 20) 

Thus, p^^' is independent of y and z. Now (jA 191) with boundary condition derived from 
(|A8|) may be readily integrated, 

y(o) ^ _p(o)„^ + ,,^^(0) ^ _p(o)yp + u,/^i _ a0(o))_ (A 21) 

where Up is the solution of {up)yy + {up)zz — —1 that vanishes at the wall. To determine 
the unknown function p*^"^ we use (IA6[) . Integrating over the cross-section and using (jASp 
we get 

MW = -p(°)up + u,/(l-a0("))=u„ (A22) 

which determines p\ . Substituting this in (lA 21 [ ) we obtain th e axial flow velocity which 
is identical to the lubrication theory solution ([Ghosall 12002 c|) applied to the case of an 
axially varying slip velocity m,/(1 — a^^*^^) along the channel wall: 



u(") 



u* 



1 - a</>(0) 



Up{y,z) 



(A 23) 



Since ^(°) = (j>''°\ (jXTt and (|XT3| show that the two dimensional vector {Y''°\Z^°'>) 
is both solenoidal and irrotational in the y — z plane, and, on account of the zero flux 
condition, must within the domain. Therefore, 

f(o)=zW=0. (A 24) 

A. 2. First Order 
From (|X2l) 

+ [(^(0) + y(o))</>(o)]^ + [(^(0) + Z(0))^(")],. (A25) 

Integrating (JA25I) over the cross-section of the capillary and using the boundary condi- 
tions (jA8p and (|A10I) we get a solvability condition for (JA25I) 



0V°^-H[(w(o)+X(o))0(o)]4 = O (A26) 

which reduces to 

4°) + [K + (1 - a^^°Y^W% = (A 27) 

on using (|A16p and (jA23[) . On using the solvability equation (|A27p in (JA25P we get 

Pe~'H'J + <Pi'J] = [("'°^ + X^"^ - "^"^ - X("^)^^% - uf ^(«) (A28) 

where we have used the continuity equation (|A6|) to replace Vy + Wz with ~ui . 
Equation (jA 28p is readily integrated. If the constant of integration is expressed in terms 
of (j)^^^ we get 

</,(!) = ^(1) -aPe u,G4>'-''m°^/[up (1 - a^'")) ] (A 29) 
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where G is defined by the equation Gyy + Gzz — Up — Up with the condition that G has 
zero flux at the walls and is normalized so that G = 0. From (lA 1[) 

[(1 - a(/.(°))y(i)]j, + [(1 - a0("))Z(i)]^ = a[(/.(i)x(o)]5 - [(1 - a<j>^'^^)X^% (A30) 

where we have used (|A24p . Integrating over the cross-section and using the boundary 
condition gives the solvability condition 

ai^'^'^X'^% - [(1 - a^(°))X(i)]^ = 0, (A31) 

which may be integrated with the far field boundary conditions to give 



(l-a(/)(0))2 
A. 3. Second Order 



(A 32) 



From 



Integrating over the cross-section and using the boundary conditions we derive 



+ 0V + [u^°^^^^ + X(o)(/.(i) + X(i) 



= Pe" 



''W 



(A 33) 

(A 34) 



A. 4. Reconstitution 

On combining (jA34[) with (jA27[) after substituting for the lower order terms from (|A 16[) . 
(|X23| and (|X32l) we get 



e^^i^) + e 



2^(?) 



1 



-eu* 



^(1) 



(1_Q,0(O)) 

^(0)0(1) 



^kPe 



■(l-a^(o))2 



'l-a0(o) (l-a(/)(o))2 



e2pe-i0g)=O, 



(A 35) 



which may be put in an alternate form that is asymptotically equivalent to it for all 
terms up to order e^: 

a2;. 

+ O{e^)^0. (A 36) 



1 



u* 



1 — aq 



a kPe u^- 



'(l-a(/ 



Pe' 



Here k — —Gup/uz, is a constant that depends solely on the geometry of the channel 
cross-section. On transforming back to dimensional variables p.23p follows. 



